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We  consider  here  the  problem  of  calculating  the  Q  of  a  very  overmoded,  (wavelength  <  <  cavity 
size)  irregularly  shaped  resonant  cavity  due  to  absorption  of  the  electromagnetic  radiation  at  the  walls. 
We  assume  further  that  the  volume  inside  the  cavity  can  be  either  vacuum,  or  else  partially  filled  with 
anisotropic  inhomogeneous  dielectric  or  plasma.  Both  the  scale  length  of  the  dielectric  and  the  radius 
of  curvature  of  the  walls  are  assumed  much  larger  than  the  radiation  wavelength.  However,  near  the 
cavity  walls,  a  vacuum  is  assumed.  Our  main  interest  here  is  in  applications  concerning  magnetically 
confined  plasmas.  For  instance  consider  a  tokamak^  containing  a  hot  plasma  which  radiates  at  the  cyclo¬ 
tron  frequency  and  its  harmonics.  An  important  issue  is  how  much  of  this  cyclotron  radiation  is 
absorbed  by  the  walls  and  how  much  is  reabsorbed  by  the  plasma.  Generally  this  can  be  calculated  with 
a  ray  tracing  code  where  many  rays  are  followed  and  the  absorbtion  and  emission  are  calculated  along 
each  ray  path.  However  if  the  plasma  is  optically  thin,  as  it  would  be  at  the  higher  harmonics,  and  if 
the  reflection  coefficient  at  the  wall  is  near  unity  the  rays  would  have  to  be  followed  for  long  dis’ances 
before  one  could  see  how  much  energy  is  deposited  in  the  wall.  ^ 


v  ^  K 

In  this  paper  we  utilize  the  ergodic  theorem  to  calculate  the  wall  absorbtion  for  the  case  where  the 
wave  makes  many  bounces  and  before  it  is  absorbed.  If  the  local  dispersion  relation  is 

«  -  (1) 

where  the  F  is  real,  the  ray  equations 


dr  dF  d  k  _  &F 

dt  "  dk  ’  dt  d_r  Ul 

are  relevant.  It  is  evident  from  Eqs.  (2)  that  the  ray  equations  are  a  Hamiltonian  system  for  which 

Fijs.r)  plays  the  role  of  the  Hamiltonian,  and  Oc,_r)  are  momenta  and  coordinate  variables.  Thus 
Eqs.  (2)  are  subject  to  the  same  phenomena  of  ergodicity  onset  and  ergodic  motions  of  other  Hamil¬ 
tonian  systems.  Here  we  consider  the  case  where  the  solution  to  Eqs.  (2)  is  ergodic1-2  with  F  the  only 
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constant  of  the  motion  of  the  ray  equations.  (We  note  that  even  in  plasma  devices,  such  as  tokamaks, 
where  the  plasma  has  a  high  degree  of  toroidal  symmetry,  the  conducting  walls  surrounding  the  plasma 
commonly  have  corrugations  and  are  partly  composed  of  baffels  and  limiters.  Thus  this  ergodicity 
assumption  is  probably  well-justified  for  waves  which  experience  reflections  from  the  walls.)  In  the  case 
of  ergodic  ray  motions  it  is  expected  that  the  time  average  over  the  motion  of  a  wavepacket  can  be 
obtained  by  a  phase  space  average 

lim  t-1  f  (/), _r(f)]<*  -  f  dj k  f  drf(k,_r)  g(k  ,_r)  (3) 

T  — **  U  - 

where  g  is  any  function  of  k  and  r  and  /  is  the  microcanonical  distribution 


f(k  ,_r)  —  8(w 


dk  drb  (o>  -  F(.k,r)) 


(4) 


On  the  basis  of  (4),  it  is  expected  that  for  an  eigenmode  with  resonant  frequency  «■>  the  wave  energy 
density  in  k,  _r  is  given  by3 

-  Iff  w(kt  r)  dkdj\  -■ rr --- 

1  Sfdkdr 


~  FUf.r » 


(5) 


8 (to  -  F(k,_i )) 

if  the  solution  of  the  ray  equations  is  ergodic.  This  is  analogous  to  the  microcanonical  ensemble  in  sta¬ 
tistical  mechanics  where  all  phase  space  points  on  the  energy  surface  (analogous  to  the  constant 
F(k._r)  surface)  are  equally  likely.  In  writing  Eqs.  (4)  and  (5)  we  have  in  mind  a  cavity  filled  with  a 
magnetized  plasma  or  anisotropic  dielectric.  Also  we  consider  only  electromagnetic  waves  which  pro¬ 
pagate  freely  from  the  vacuum  into  the  plasma  and  visa  versa.  In  this  case  the  two  polarizations  of 
plane  electromagnetic  waves  propagating  in  a  given  direction,  with  given  frequency  are,  in  general, 
non-degenerate.  That  is  they  have  different  values  of  i  k  I.  Thus  the  dispersion  relation  a>  *■  FM.j;) 
will  possess  these  two  solutions.  Hence  in  this  case,  Eq.  (5)  automatically  gives  the  distribution  of 
energy  in  the  two  polarizations.  In  the  case  of  an  unmagnetized  plasma,  isotropic  dielectric  or  vacuum, 
the  two  independent  polarizations  are  degenerate,  that  is,  they  have  the  same  value  of  \k  |.  In  this 
case  we  supplement  Eq,  (5)  with  the  information  that  both  polarizations  are  equally  likely.  One  expects 
that  Eq.  (5)  holds  in  the  limit  of  k  —  <»■,  however,  it  is  not  yet  known  how  good  Eq.  (5)  is  for  finite  k. 
In  this  respect,  recent  numerical  experiments  are  of  interest4  (typically  k  L  —  10  in  these  experiments 
where  L  is  a  scale  length).  We  shall  be  interested  in  applying  Eq.  (5)  in  cases  where  the  wave  experi- 
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ences  absorption  at  the  conducting  boundaries  upon  reflection.  Under  these  circumstances,  we  require 
the  following  two  conditions  for  the  validity  of  Eq.  (5).  Condition  A:  The  wall  absorption  in  following  a 
typical  ray  must  be  slow  enough  that  the  ray  wanders  over  a  representative  region  of  the  “ 

constant  surface  before  significant  absorption  takes  place.  Condition  B:  The  dispersion  relation  E(k,j) 
must  not  exhibit  any  resonances  (i.e.,  solutions  with  k  —•  °°).  Condition  B  is  necessary  in  order  that 
the  integral  in  the  denominator  of  Eq.  (4)  be  finite. 


We  now  calculate  the  damping  o:'  an  ergodic  eigenmode  due  to  finite  wall  conductivity.  The  Q  of 
the  cavity  is  defined  as 


<?- 

where  tj,  is  the  surface  skin  resistivity  tj 


-//  dk  dr 

$V,\J,\2d2r 

,  —  <j*0u/2 <r)1/2,  cr  is  the  wall  conductivity  which  may  be  a 


function  of  position  on  the  wall,  Jx  is  the  rms 


surface  current  density  vector,  and  the  integral,  d2r. 


is  taken  over  the  boundary  of  the  cavity.  Here  we  wish  to  calculate  \  JS  1/  (and  hence  Q). 


We  recall  that,  by  assumption,  near  the  cavity  walls  the  plasma  density  is  zero.  In  this  region 
Fijt.j)  "  I  k  |c,  where  c  is  the  speed  of  light,  and  from  Eq.  (5)  W(k,z)  —  -  |fc|c].  Thus,  near 

the  walls,  H'Qc,  r)  is  independent  of_r  and  the  direction  of  k.  To  treat  this  case  we  recall  that  k  «  L 
and  consider  a  plane  wave  in  vacuum  which  is  incident  on  a  plane  conductor.  We  also  assume,  as  men¬ 
tioned  previously,  that  in  the  vacuum  region  near  the  wall  W{k,r)  is  the  same  for  each  v.'  the  two 
independent  polarizations.  Let  6  be  the  angle  of  incidence,  and  <t>  the  angle  that  the  magnetic  field  vec¬ 
tor  H  makes  with  the  plane  of  incidence.  Then,  from  Maxwell’s  equations  and  the  condition  that  E 
tangential  vanish,  the  current  7,  created  by  this  plane  wave  is 

7,-  2//cos*a  +2//sin*cos#£  (7) 

where  a  and  £  are  unit  vectors  in  the  surface  of  the  conductor,  respectively  perpendicular  to  and  in  the 

plane  of  incidence.  Averaging  over  <t>  and  the  solid  angle  of  the  incident  waves  we  have  that 

\JS  I2 1  *  J*0  f  (4//2(cos2$  +  sinfycos2#)]  (8) 
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where  d(l  —  2ir  sirtJ  d9,  the  integral  over  dCl  is  from  0  —  0  to  9  —  n/2  (i.e.,  only  over  incident 
waves),  and  H 2  may  be  related  to  the  average  energy  density  of  the  eigenmode  in  the  vacuum  near  the 
wall 


1  ,  rr  ,  x  k  f  siw  -  l/c  Id  dk 

**  f  f  W{k.r)  dk  dr  - . 

2  Jj  8  [w  -  /r(fclr)]  dj(dr 


(9) 


where  V  denotes  the  volume  enclosed  by  the  cavity.  From  Eq.  (8),  Us  |J  -  4/3  H J.  Thus 


Q 


3mo<’3 

16u  w  £)  Tjjt/2rK 


8 [tu  -  F(fc,r)]  d/edr. 


(10) 


For  vacuum,  F(k,r)  —  |*|c  and  Eq.  (10)  yields  Q  —  {3^owl  isd2r]  1  where  the  two  indepen¬ 
dent  polarizations  have  been  taken  into  account. 


Similar  arguments  to  those  above  can  be  applied  to  Helmholtz  equation  in  two  dimensions  (x  and 
_y),V2^  +  lcty  —  0  with  the  eigenfunction  zero  on  the  boundary  with  the  result  that 

<  (80/cb))2  >  —  k2f$2dxdy  /  f  dxdy,  (11) 

where  d/fb)  denotes  the  normal  derivative,  and  <  . . .  >  denotes  the  average  over  the  boundary.  This 
expression  has  been  checked  in  numerical  experiments  by  McDonald  and  Manheimer5  using  the  model 
of  Ref.  6.  Reasonable  agreement  was  obtained,  thus  giving  an  indication  of  the  correctness  of  our 
approach.  To  obtain  (11)  consider  a  plane  wave  of  amplitude  incident  on  a  planar  boundary  where 
—  0.  This  leads  to  a  normal  derivative  of  at  the  boundary  of  amplitude  2$  k  cos 9,  where  9  is  the 
angle  of  incidence.  Averaging  (2 ii>  k  cos 9)7  over  0  for  incident  waves  (0  <  0  <  n/2)  then  yields  Eq. 
(11). 


If  the  wave  has  non-zero  damping  rate  ■y(k,_r)  in  the  plasma  then  the  cavity  Q  is  (0ht'  +  C?,-1)-1 
where  Qf  is  the  Q  resulting  from  absorption  in  the  plasma  and  Qw  is  given  by  Eq.  (10).  Thus  the  frac¬ 
tion  of  energy  absorbed  by  the  wall  is  QP(Q,  +  Q*)-'.  The  Q  resulting  from  the  plasma  losses  is 


Q, 


(12) 
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where2 

_  JfdiLte  >(&..£)  ®(“  ~  F(Ic,z)) 

y  “  - -jc-i i - - - .  (13) 

J  J  die  4r  *(*>-  /’'(&, _r)) 

This  may  be  useful  in  calculating  the  amount  of  energy  emitted  by  cyclotron  radiation  which  is 
absorbed  by  the  wall  in  a  thermonuclear  fusion  device. 
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